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Abstract 

Relationships between inter-cluster synchronization phenomena and external noise are studied on the basis of noise level-free 
analysis. We consider a mean-field model of ensembles of coupled limit cycle oscillators with two natural frequencies, which are 
subjected to external white Gaussian noise. Using a nonlinear Fokker-Planck equation approach, we analytically derive the order 
parameter equations associated with nonequilibrium phase transitions in the thermodynamic limit. Solving them numerically, we 
systematically investigate the model parameter dependences of the appearance and disappearance of synchronization phenomena. 
Demonstrating bifurcations from chaotic attractors in the deterministic limit to limit cycle attractors with increasing noise intensity, 
we confirm the occurrence of nonequilibrium phase transitions including inter-cluster synchronization induced by external white 
Gaussian noise. 

Keywords: Noise induced synchronization, Mean-field model, Nonlinear Fokker-Planck equation, Nonequilibrium phase 
transitions, Stochastic limit cycle oscillators 
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1. Introduction 

Effects of noise on synchronization phenomena in oscilla- 
tory systems have recently attracted much attention from many 
researchers. Synchronization phenomena are ubiquitous ones 
observed in various fields of natural sciences QJJ] . In neuro- 
sciences, neurons in the basal ganglia exhibit more synchrony 
in Parkinson's disease than in normal state, suggesting that neu- 
ral information coding is closely related to the synchronization 
phenomena Jj]- For such a reason, to study how the noise ex- 
erts its influence on the structure of synchronization will be of 
paramount importance from the viewpoint of nonlinear dynam- 
ical control involving changes in synchrony of oscillatory sys- 
tems. 

On one hand, someone would intuitively suppose that the 
presence of noise might deteriorate the degree of synchro- 
nization of oscillatory systems. Breakdown of synchroniza- 
tion due to external noise of an ensemble of limit cycle os- 
cillators was reported JsJ]. On the other hand, the opposite 
phenomena of noise-induced synchronization are becoming 
an active field of the study of nonlinear dynamical systems 
flSHSHS- Noise induced synchronization in coupled ex- 
citable systems/active rotators is investigated both analytically 
and numerically JfJ. Among analytical studies on syn- 
chronization phenomena of ensemble of limit cycle oscillators 
is the phase reduction analysis. Such type of studies revealed 
the effects of common noise on synchronization of uncoupled 
oscillators [7] and uncoupled two populations of oscillators J^l. 
These studies, however, are restricted to the case with weak 
noise. Numerical approaches to noise induced synchronization 
of coupled limit cycle oscillators might include the possibili- 
ties that what are so obtained are subjected to a finite size effect 



without taking the thermodynamic limit as in Ref . |0] . 

To avoid such limitations, one may take nonlinear Fokker- 
Planck equation approaches, which are very closely related to 
the study of (thermal) equilibrium phase transitions. Present- 
ing the validity of an H theorem for a nonlinear Fokker-Planck 
equation for a stochastic system of mean-field coupled over- 
damped oscillators, Shiino studied statistical behaviors of the 
system exhibiting equilibrium phase transitions [10] (see also 
Ref. Ill U ll2fl ). Furthermore, in the case of nonequilibrium 



phase transitions, one can take advantage of using nonlinear 
Fokker-Planck equations to exactly derive the time evolution of 
the order-parameter equations in the thermodynamic limit for 
systems of nonlinearly mean-field coupled oscillators. There, 
the nonlinear dynamic aspects of phase transitions were investi- 
gated, yielding the occurrence of chaos-nonchaos phase transi- 
tions and those including intra-cluster synchronization induced 
by external noise 01311 . Purely noise induced phase transitions 
involving chaos-nonchaos phase transitions in an ensemble of 
limit cycle oscillators were also explored II 1 411 . A nonlinearly 
coupled system with time delay was found to exhibit various 
types of nonequilibrium phase transitions lfl5ll . 

In the present study, we apply the nonlinear Fokker-Planck 
approach to noise induced inter-cluster synchronization phe- 
nomena of coupled limit cycle oscillators. Dealing with a solv- 
able model based on the mean-field concept to derive order pa- 
rameter equations II 1 3L 1 14 , we study the effects of noise 
on synchronization in the thermodynamic limit and nonequilib- 
rium stationary states, for which numerical approaches to solv- 
ing stochastic differential equations do not properly work. Our 
results show the appearance of nonequilibrium phase transitions 
involving inter-cluster synchronization. Since the model of os- 
cillators is based on that of analog neural networks, behaviors 



Preprint submitted to Physica E 



January 4, 2011 



of noise-induced synchronization in the system will be of prac- 
tical interest. 



2. Model and nonlinear Fokker-Planck equation approach 

We consider a system of coupled limit cycle oscillators con- 
sisting of two subsystems. The oscillators are coupled via non- 
linear global interactions and subjected to independent addi- 
tive and multiplicative external noise. The model dynamics we 
are concerned with is described by a set of Langevin equations 
Q: 
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where zjv (fi — x,y) (a — 1, 2 (3 = 1)) are the 2D-oscillators at 



site i, da , b^ ,V) are constants and F^'(-) are bounded functions 
specifying the nonlinear couplings. The mean-field coupling 
strengths J^fXt) are given by 
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In cluster 1 and 2, we postulate the different constraint parame- 
ters as = a i + w i tn <^ being essentially responsible 
for the difference of natural frequencies between the oscillators 
in the two subsystems. The inter-cluster coupling strengths are 
controlled by e. This parameter takes any real constant value, 
without constraining to weak connections. The external noise 
rj^{t), g£pij(f) are of white Gaussian type, 
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Note that not common but independent noise is introduced and 
rj^l (f), Sjy(0 are a l so independent each other. 

In the thermodynamic limit n a — > 00, influence with Eq. <£3j 
on any one oscillator from the others almost surely converges 
according to the law of large numbers, satisfying the validity of 
the self-averaging property. Introducing empirical probability 



density P(t,z\ , Z{ , Z^ ,Z% ), the mean-field coupling terms in 
Eqs. ([TJ and © are expressed in terms of 

xP(t,z ( f,z^\z ( f,zf), (4) 
{F { f) = j dzfdzfdzfAzfFfi^z^+b^zf) 

xP(t,z ( f,z^\z 2 x \z^). (5) 

The total number of dynamic variables of the system is con- 
sequently reduced from 2(«i + ni) to 4 in the thermodynamic 
limit. Then one has a set of the Langevin equations, 
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where the noise is subjected to white Gaussian noise, 

<^(0> = 0, 
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The nonlinear Fokker-Planck equation lfl6ll corresponding to 
Eqs. © reads 
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A Gaussian probability density satisfies Eq. © as a special 
solution. Since the H theorem 111 311 ensures that the probability 
density satisfying Eq. © converges to the Gaussian-form for 
sufficiently large times, we represent the Gaussian probability 
density as 



P G (t,z ( ?,zf\z ( ?,z ( ?) = 



(2n) 2 VdetC G (f) 



,(*) 



exp 



,O0\ 



-~4qV>c 



,(*) 



<Z? ) )c,Z? ) -<Z? ) ) ) 



_ , (x) (>•) (x) 

C Gij (t) = (s t Sj) G , 

where (-) G denotes expectation over P G . Now that Eqs. © and 
© are described in terms of the first and second moments, one 
obtains a set of closed ordinary differential equations involving 
at most second moments. The time evolution of each moment 
is calculated from Eq. (0 and we write the moment equations 
as 

d<zr>c 



= -af{zf) + JZ(K')« + eCi.<Ci >° ( 8 ) 



7<M) 



df 
df 

d(«r4 V) ) G 

d7 



-2fl£\ur>o+2D<£ 



-(a^+4 y) )< M r4 V) )c 



(9) 
(10) 



where 

(a,/3,n,v) = (l,l,x,y), (2,2,x,y), (l,2,x,x), 
{\,2,y,y),{\,2,x,y),{\,2,y,x). 

Note that (u^u^} a — > (f — > °°), implying that the covariant 
components of C G take zero in the stationary states. 

For observing qualitative as well as quantitative dynamical 
behaviors of the system, we proceed to solve numerically the 
above set of equations. We specify the coupling function as 
F^\x) = sin(x), for simplicity. Then, (F^) G and (F^ ) G are 
calculated as 



x exp 
1 



2 v a „ 2 

<F^% = J - 5 cos (2b^\z?) G + 2b^\zt\) 
x e x V (-2b^(uf) G -2b^\uf ) G ). 



(ID 



(12) 



3. Nonequilibrium phase transitions and noise induced 
inter-cluster synchronization 

The dynamic system of the order parameter equations written 
by Eqs. (T8b-(fT0b together with Eqs. dTTb and ( TTZl i is nonlinear, 
allowing various types of bifurcations with changes in the pa- 
rameters. For the definition of synchronization, let us consider 
relationships between behaviors of statistical variables in Eqs. 
(l8ll-(fT0t and those of individual oscillators in the system. 

First, we address the issue of synchronization phenomena of 
the coupled oscillators in the stochastic case. On one hand, 
in deterministic systems of the coupled limit cycle oscillators, 
the synchronization phenomena are well-defined by the condi- 
tion that all of the variables of the system undergo a periodic 
solution with a common time period. On the other hand, in 
stochastic systems of a finite number of coupled limit cycle os- 
cillators with external noise, each oscillator behaves randomly 
under the influence of noise. Since the Fokker-Planck equations 
in this case are linear, the probability densities of the systems, 
in general, exhibit equilibrium ( i.e. fixed point type ) probabil- 
ity densities for sufficiently large times. This implies that the 
order parameters do not oscillate even if individual oscillators 
behave periodically in the deterministic case. It seems to be dif- 
ficult to appropriately define the synchronization phenomena in 
this situation. 

To consider the synchronization phenomena in systems of 
stochastic coupled limit cycle oscillators without the problem 
mentioned above, we have introduced the concept of taking the 
thermodynamic limit based on a mean-field model. In general, 
probability densities of the limit cycle oscillator systems might 
be multimodal under the influence of weak external noise. To 
avoid the difficulty of dealing with multimodal probability den- 
sities, our model is proposed so that once external noise is in- 
troduced to the system, the model constrains the form of the 
probability densities to the Gaussian ones for sufficiently large 
times. Note that the probability densities may be allowed to 



vary in time. In this case, the order parameters of the system 
may oscillate periodically in time and we can define the syn- 
chronization phenomena according to their behaviors. In par- 
ticular, the oscillations of the mean values have an important 
meaning from the view point of analogy with the spontaneous 
magnetizations in spin models and possibly researches in neu- 
rosciences. 

For these reasons, we consider synchronization phenomena 
of the oscillators in each cluster a, i.e., intra-cluster synchro- 
nization and between cluster 1 and 2, i.e., inter-cluster synchro- 
nization. We define intra-cluster synchronization phenomena 
as (i) the case that "all of the mean values over the oscillators 
in a cluster periodically oscillate to form limit cycle attractors, 
regardless of the behaviors of the variances," and inter-cluster 
synchronization phenomena as (ii) the case that "(i) occurs as 
a whole system." We focus our attention on the inter-cluster 
synchronization phenomena. 

For simplicity, we only treat the Langevin noise case, i.e., 
D x — D y = 0, and do not study inter-cluster chaotic synchro- 
nization in this paper. From Eqs. the variances turn out 
to be constant for sufficiently large times. To investigate the 
effects of the Langevin noise on the appearance and disappear- 
ance of the inter-cluster synchronization phenomena in the ther- 
modynamic limit and nonequilibrium stationary states, we con- 
fine ourselves to the bifurcations of limit cycle attractors with 
changes in the noise intensity. To determine whether an attrac- 
tor is of limit cycle type or not, we employ two criteria: (a) the 
largest Lyapunov exponent (LLE) of the attractor estimated nu- 
merically nearly equals zero and (b) all of the mean values have 
a common time period. Note that both of them exclude the case 
where the attractor is of ergodic torus type. 

Numerical calculations were performed with the fourth-order 
Runge-Kutta method. The model parameter values are a. — 



0.5, af = 1.0, b\ 



(x,x) 



0.5, b\ 



1.5, b 



(y,x) 



-12.0, b 



0',.v) 



-1-0, J$ = J% = 18-0, Jg = J% = 30.0, af = 0.5 + 
flj - 1.0 + d^. Nonequilibrium phase transitions involving 
the inter-cluster synchronization of the limit cycle oscillators 
are systematically investigated with changes in the inter-cluster 
coupling strength e and difference of constraint parameters <$^0 
, and the Langevin noise intensity . 

To show the appearance of the inter-cluster synchronization 
due to the effects of noise, we begin with investigating behav- 
iors of the system in the deterministic limit, especially the case 
where a chaotic attractor appears (Fig. [TJ. With increase of the 
Langevin noise intensity D (x \ the attractor in this case changes 
into limit cycle type via torus type, suggesting the appearance 
of the inter-cluster synchronization. Further increase of the 
noise intensity leads to the negative LLE, implying fixed point 
type attractors. These results are qualitatively consistent with 
those observed in Ref. [9], which numerically solved a set of 
Langevin equations to be subjected to the finite size effect. 

For clearly demonstrating the effects of noise on the inter- 
cluster synchronization, let us show the temporal evolution of 
the mean values (z\ ) a , {z\ * )g of the oscillators in cluster 1 and 
2 (Fig. [2J. It is easily seen that the strength of noise above a 
certain level causes the synchronization between the clusters, 
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Figure 1 : The largest Lyapunov exponents plotted against the Langevin noise 
intensity Br e ' . In the deterministic limit, the LLE has a positive value, imply- 
ing a chaotic attractor. With an increase of the noise intensity, the LLEs take 
nearly equal zero, and the attractors changes into limit cycle type via torus type 
(denoted by LC and T, respectively). We set Eft' = = 0. 



although both of the oscillators chaotically behave in the deter- 
ministic limit. 



4. Concluding remarks 

Using a noise level-free analysis, we have shown the relation- 
ships between inter-cluster synchronization phenomena and ex- 
ternal noise, especially the occurrence of the inter-cluster syn- 
chronization induced by the Langevin noise. Dealing with a 
mean-field model of ensembles of coupled limit cycle oscilla- 
tors with two natural frequencies under the influence of external 
white Gaussian noise, we have analytically derived the order 
parameter equations in the stationary states, taking advantage 
of the ingredients of the self-averaging property in the thermo- 
dynamic limit together with an H theorem that ensures the con- 
vergence of the probability density to the Gaussian-form for 
sufficiently large times. Solving the order parameter equations 
numerically, we have systematically investigated the nonequi- 
librium phase transitions with changes in the inter-cluster cou- 
pling strength e and difference of constraint parameters S^', and 
the Langevin noise intensity D^K The results have shown vari- 
ous interesting bifurcations including the inter-cluster synchro- 
nization and chaotic attractors induced by noise. 

Details of the relationship between the parameters involv- 
ing the inter-cluster coupling strength e as well as difference of 
constraint parameters S <p) and the behaviors of the inter-cluster 
synchronization, and furthermore the effects of the multiplica- 
tive noise on the inter-cluster synchronization will be reported 
elsewhere. 
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Figure 2: The time evolution of the mean values (z\ )a, (z 2 )a of the oscillators 
of the oscillators in the two clusters, corresponding to the changes from the 
chaotic to limit cycle attractors in Fig. \T\ While the oscillators do not behave 
as synchronous in the deterministic limit (the upper panel), an appropriate level 
of the Langevin noise intensity induces the inter-cluster synchronization (the 
lower panel). 
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